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Abstract: The CALICE collaboration is studying the design of high performance electtomag- 
netic and hadronic calorimeters for future International Linear Collider detectors. For the hadronic 
calorimeter, one option is a highly granular sampling calorimeter with steel as absorber and scin- 
tillator layers as active material. High granularity is obtained by segmenting the scintillator into 
small tiles individually read out via silicon photo-multipliers (SiPM). A prototype has been built, 
consisting of thirty-eight sensitive layers, segmented into about eight thousand channels. In 2007 
the prototype was exposed to positrons and hadrons using the CERN SPS beam, covering a wide 
range of beam energies and angles of incidence. The challenge of cell equalization and caUbration 
of such a large number of chaimels is best vaUdated using electromagnetic processes. The response 
of the prototype steel-scintillator calorimeter, including linearity and uniformity, to electrons is 
investigated and described. 

Keywords: Calorimeter; electromagnetic shower; Silicon Photomultiplier. 
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1. Introduction 



A new generation of calorimeters that exploit unprecedented high granularity to reach excellent jet 
energy resolution is one of the main R&D goals towards the future International Linear Collider 
(ILC) |jl|]. The particle flow (PFLOW [|2[ ||, ^) algorithm favors single particle separability over 
single particle energy resolution in the attempt to improve the overall jet energy resolution. Typical 
single hadronic showers in the 10-100 GeV range are best separated in a hadronic calorimeter with 
cell size of the order of 3 x 3 cm^ [Q]. In addition, fine longitudinal segmentation is required for 
PFLOW algorithms to be effective. 

The CALICE collaboration |j5|] is studying several calorimeter designs for experiments at the 
ILC. With the first generation of prototype detectors new readout technologies have been estab- 
hshed for highly granular calorimeters and the stability of these detectors has been demonstrated. 
Furthermore, a unique set of data has been collected to study hadronic showers at low and medium 
energies in detail with high resolution longitudinal and transverse sampling. 

This paper focuses on the prototype of an analog hadron calorimeter (AHCAL) consisting 
of 38 layers of highly-segmented scintillator plates sandwiched between 2 cm thick steel plates. 
Each scintillator tile is an individual calorimeter cell read out by a silicon photo-multiplier (SiPM). 
Details on the calorimeter structure, calibration and readout electronics are given in Section |2[ 

Tests using particle beams have been conducted in order to evaluate the performance of the 
highly granular calorimeters built by CALICE. 

In 2007, the whole detector with 38 active layers was commissioned and exposed to muon, 
positron and pion beams in the energy range 6 GeV to 80 GeV provided by the CERN Super Proton 
Synchrotron (SPS \^), on the H6 beam hne. In 2008, the AHCAL together with the ECAL and 
TCMC were moved to the FNAL Meson Test Beam Facility (MTBF [|l2|]) to take data in the 1- 
6 GeV energy range over the course of two years. 

Ongoing data analyses will quantify the energy and spatial resolutions of the prototype for 
hadrons, and will continue the validation and further development of existing models of hadronic 
showers, e. g. the various Geant4 [|1^] physics lists. They will also be important for the experi- 
mental validation of the PFLOW approach [|l^]. The studies in this article focus on the calibration 
and performance of the device when exposed to electrons and positrons. 

In Section ^, the AHCAL is described, and in Section ^ the electromagnetic calibration proce- 
dure is discussed. The CERN test beam experiment is described in Section |^ Results on calorime- 
ter response to positrons are given in Section ^, followed by uniformity studies in Section ^ Con- 
clusions are reported in Section 0. 



2. Prototype calorimeter 

The AHCAL is a sampling calorimeter with alternating 2 cm thick steel plates and highly-segmented 
scintillator-based active layers. The single calorimeter cell is a scintillator tile read out via a SiPM. 
The scintillator tiles are 0.5 cm thick and have a size of 3 cm x 3 cm in the 30 cmx 30 cm core region 
and increase to 6 cmx 6 cm and 12 cmx 12 cm in the rings surrounding the core. A sketch of one 
AHCAL module, as well as a picture of an open module showing the arrangement of scintillator 
tiles, are shown in Figure [|. A wavelength shifting fiber is embedded in the tile, which collects the 
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AHCAL MODULE 




Figure 1. Sketch of one AHCAL module (left). The scintillator tiles with SiPM readout are embedded in a 
steel cassette. The SiPM signal is routed to the VFE electronics located on the right side. The calibration and 
monitoring board (CMB), located on the left side, provides UV LED light for calibration and the CAN-BUS 
readout for temperature sensors located inside the cassette and on the electronics (red dots). Picture of one 
active layer of the CALICE AHCAL prototype (right). 

scintillation light and guides it to the SiPM. The other fiber end is pressed against a 3M reflector 
foil. The details of one core tile are shown more clearly later in this paper, in the left side of Fig- 
ure l^. The four sides of each tile are matted by a chemical treatment providing a white surface 
that serves as a diffuse reflector. The two large faces of the tile are not individually coated, instead 
a large 3M reflector foil is glued to each side of the metal cassette hosting all tiles providing reflec- 
tivity via an air contact. The imperfect reflective coating of the tile edges is responsible for about 
2.5 % light cross-talk between neighboring cells of 3 cm edge size. 

SiPM devices from the MEPhl/PULSAR group have been used, which have an active area of 
LI mmx LI mm containing 1156 pixels, each 32 /xmx32 /im in size. SiPMs are operated with 
a reverse bias voltage of ~50 V, which lies a few volts above the breakdown voltage, resulting 
in a gain of ~10^. A poly-silicon quenching resistor on each pixel is used to quench the Geiger 
discharge. The resistor values vary between 0.5 Mil and 20 Mil for the various batches of SiPM 
produced. Larger resistor values have been favored as they yield a longer pixel recovery time up to 
1 jUs. In this way a pixel cannot be fired multiple times during one scintillator light pulse, though 
making it easier to monitor the SiPM response curve with LED light. More details on the SiPM 
working principle and its properties are given in [§, 0]. 

The active layers are referred to as modules, and the sum of active and passive material adds 
up to a total depth of 5.3 nuclear interaction lengths (A,)- A more detailed description of the 
AHCAL prototype structure is given in [j8|]. The analog SiPM signal is routed to the very-front-end 
(VFE) electronics where a dedicated ASIC chip is used for multiplexed readout of 18 SiPMs. 
The integrated components of the ASIC chip allow to select one of sixteen fixed preamplifier gain 
factors from 1 to 100 mV/pC, and one of sixteen CR-RC^ shapers with peaking times from 40 to 
180 ns. 
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Since the AHCAL was the first detector to employ such a large number of SiPMs, a specialized 
system for monitoring the long-term stability and performance of the photodetectors was required. 
In order to monitor the SiPM response function in-situ, a versatile UV LED light distribution 
system was developed [[l^]. A calibration and monitoring board (CMB) connected to each module 
distributes UV light from an LED to each tile via clear fibers. The LEDs are pulsed with 10 ns 
wide signals steerable in amplitude. By varying the voltage, the LED intensity covers the full 
dynamic range from zero to saturation (about 70 times the signal of a minimum-ionizing particles). 
Furthermore, the LED system monitors variations of SiPM gain and signal response, both sensitive 
to temperature and voltage fluctuations. The LED light itself is monitored with a PIN photo-diode 
to correct for fluctuations in the LED light intensity. 

3. Calibration procedure 

One of the aims of the tests discussed here is to establish a reliable and robust calibration chain. 
This requires measurements with beam particles and with light from the LED monitoring system. 
The calibration chain is summarized in the following steps: 

• calibration of the cell response and cell-to-cell equalization; 

• monitoring of the SiPM gain and corrections for the non-linear response; 

• calibration to an energy scale (in GeV) with electromagnetic showers. 

The cell-by-cell calibration, and with that the equalization of all cell responses, is performed 
using minimum-ionizing particles (MIPs) provided by a broad muon beam with an approximately 
Gaussian profile with a width of about 30-40 cm, illuminating all cells in the detector. For each 
cell, a calibration factor, Cf^^^, is determined from the most probable value of the measured energy 
spectrum for muons in ADC units, which is extracted with a fit using a Landau function convolved 
with a Gaussian. This fit accounts for the distribution of energy loss of muons in the scintillator 
tiles as well as for contributions from photon counting statistics and electronic noise. The combined 
systematic and statistical uncertainty for these fits was typically on the order of 2%. The muons are 
generally parallel to the beamline and perpendicular to the detector front face. In this way all cells 
can be calibrated at the same time, minimizing the impact of temperature induced variations. 

The SiPM gain and photo-detection efficiency are temperature dependent. The product of the 
two determines the SiPM response, which typically decreases by 3.7%/K. A procedure has been 
developed to correct temperature-induced variations in the calorimeter response using temperature 
measurements in each module. This procedure and its stability will be described in more detail 
in [|^]. For the analysis presented in this paper, data samples have been selected to cover a tem- 
perature range of less than 0.5 K to reduce the impact of such corrections. To account for the 
included temperature variations, the visible energy of each data set is scaled by -3.7%/K to the 
average temperature of the muon data used for calibration. 

The number of SiPM pixels, A;[pix], firing for a single cell / is related to the ADC value for the 
cefl. A; [ADC], and the corresponding SiPM gain, Cf" [ADC] by A, [pix] = A,- [ADC] /Cf" [ADC] . 
The procedure to obtain the gain of each individual SiPM is discussed in Section [3.1| . 
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The limited number of SiPM pixels leads to a non-linear response for large signals. These 
effects are corrected for by a function, /sat(A[pix]), depending on the number of fired pixels A [pix] . 
This procedure is discussed in detail in Section p.2[ 

Finally, a common calibration factor, w, scales the visible energy of electrons in each cell 
in units of MIP to the total deposited energy in units of GeV. This factor is determined to be 
w = (42.3 ± 0.4) MIP/GeV, as discussed in Section g 

Therefore, in summary, the reconstructed energy of electromagnetic showers in the calorimeter 
is expressed as 

^--[^^^] = H.[MIP/GeV] ' ^^-'^ 

where the energy of one single cell with index / is The energy, Ei, given in units of MIP is 
calculated according to 



Ei [MIP] = -lL^.f,,,{At [pix]). (3.2) 



3.1 SiPM gain and electronics inter-calibration factors 

The gain of each individual SiPM is extracted from single photoelectron spectra taken in dedicated 
runs with low LED light intensity. LED light is necessary as the best determination of the gain 
requires a single photoelectron spectrum with a Poisson mean of about 1.5 p.e. and the mean 
obtained from dark noise events is below 0.5 p.e. 

The SiPM gain, , is the distance between two consecutive peaks in the single photoelec- 
tron spectrum. A typical gain spectrum is shown in Figure ^ A multi-Gaussian fit is performed 



to the single photoelectron peaks to determine their average relative distance []17[]. The mean of 
each Gaussian function in the multi-Gaussian sum is left as a free parameter. Before fitting, a peak 
finder routine is used to set each peak mean value to the approximate location of the corresponding 
photoelectron peak. It has been seen that fixing the distance between peaks to one common pa- 
rameter reduced significantly the number of converged fits. The width of each Gaussian function 
is dominated by electronic noise, but for large number of pixels fired the statistical contribution 
becomes visible, which lead to an increase of the peak width. Accordingly, the width of each peak 
is left as a free parameter. Finally, the SiPM gain is defined from the fit result as the distance 
between pedestal (peak zero) and second peak divided by two. An additional consistency check 
is performed to ensure the distance between pedestal and first peak agrees with the defined gain 
value at better than 2%. Gaussian fits to peaks higher than the 3rd one are not directly used in 
determining the gain, but their proper description helps to improve the stablity of the fits and to 
avoid bias on the peak position. The uncertainty on the gain determination is mainly due to the fit 
and is about 2 % for fits which pass quality criteria. 

SiPM gain measurements were repeated approximately every eight hours during test beam 
operation. The SiPM gain varies with temperature and the gain measurements can be used to 
stabilize the calorimeter response over time. The temperature dependence of the SiPM gain is 



further discussed in []15p. The efficiency of the gain extraction is defined as the number of successful 
fits in one gain run divided by the number of channels which can be calibrated. About 2 % of all 
SiPMs are considered inactive because of initially bad soldering or subsequent broken connections 
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Figure 2. Single photoelectron peak spectrum taken with a SiPM in the AHCAL detector. 



to the SiPM leads. Additionally, about 0.11 % of all channels are connected to a broken LED. All 
these channels are not accounted for in the total of channels that can be calibrated. 

The efficiency of the gain extraction with one measurement run is indicative of the quality of 
the LED monitoring system, namely the small spread of LED light intensity. Figure ^ shows the 
efficiency of the gain extraction for a series of runs taken in the first three months of data taking at 
CERN and in the first three months at FNAL. Initial problems during the system commissioning 
phase led to low efficiency, but after commissioning a gain extraction efficiency of about 95 % per 
run has been achieved. The gain efficiency was also stable after transportation and throughout the 
FNAL runs. Combination of several gain runs yields calibration of more than 99 % of all cells. The 
remaining 1 % of cells are calibrated with the average of the module to which they belong. 

The measurement of SiPM gain is performed with a special mode of the readout chip, with a 
choice of high pre-amplification gain and short peaking time of 50 ns which improves the signal to 
noise ratio at the single pixel level. In contrast, the muon calibration and the physics data taking 
are performed with approximately ten times smaller electronic amplification, to optimally fit the 
available dynamic range, and about 180 ns peaking time to provide sufficient latency for the beam 
trigger. The inter-calibration factor, /,, of the chip gain between the calibration mode (CM) and 
the physics data mode (PM) along with the SiPM gain are used to determine the overall SiPM 



calibration factor, Cj [ADC], used in Eq. 3.2: 



Cf' = Gf'[ADC(CM)]//,-, 



(3.3) 



The extraction of the inter-calibration coefficients depends on the linear response of the chip 
in both modes for an overlapping range of input signals. The input signal is provided by the LED 
system injecting light into the tiles. The amplitude of the signal is varied within the linear range 
by varying the LED light intensity. The response in each readout mode is fit with a line, and 
the ratio between the two slopes is the inter-calibration coefficient for one given readout channel. 
Ideally, this factor should be a simple constant between the two chip readout modes, but it turns 
out to depend on the SiPM signal form due to the different shaping times in the two modes. For 
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Figure 3. Gain calibration efficiency (left) and electronics inter-calibration efficiency (right) over the AH- 
CAL data taking period at CERN in 2007 (red dots) and at FNAL in 2008 (open blue triangles). More than 
85.0 % of the channels could be monitored for gain and inter-calibration variation during these periods. 

longer SiPM signals (larger quenching resistor) the inter-calibration is bigger than for shorter SiPM 
signals (smaller quenching resistor). The inter-calibration factors between the chip readout modes 
range between 4 and 13. 

As with the gain, the inter-calibration extraction efficiency is influenced by the quality of the 
LED light distribution system. The inter-calibration coefficient extraction efficiencies during the 
2007 and the 2008 data taking periods are plotted in Figure ^ (right). After commissioning was 
completed, all channels with the exception of the 2 % inactive channels and the channels connected 
to a broken LED, could be inter-calibrated. For the missing inter-cahbration values the average of 
the module to which a SiPM belongs is used instead. 

The uncertainty on the inter-calibration coefficient has been estimated from the comparison of 
several runs and is found to be better than 1 %. Temperature and voltage changes do not affect this 
coefficient since it is mainly driven by the stability of the components of the readout chip and of 
the SiPM quenching resistor, all of which are stable in a range of 5-10 degrees. 

3.2 SiPM non-linearity 

Due to the limited number of pixels and the finite pixel recovery time, the SiPM is an intrinsically 
non-linear device. The SiPMs used in the AHCAL have a total of 1 156 pixels with a recovery time 
between 25 ns and 1 /is, depending on the value of the quenching resistor. 

The response function of a SiPM correlates the observed number of pixels fired, Npa, to the 
effective number of photoelectrons generated, A'pe, including cross-talk and after-pulses. The re- 
sponse of a SiPM can be approximated by the function 
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Figure 4. The SiPM non-linearity correction function, /sat- The points are the tabulated data of A^pix versus 
A^pe for one SiPM in the AHCAL. The solid line is obtained from the double exponential fit to the points, 
and it is used as a correction in this analysis. The dashed Une is the single exponential approximation from 



equation 3.5 



with A^tot the maximum number of fired pixels. This formula is a useful approximation for the case 
of uniform light distribution over the pixels and short light pulses. 

In the above approximation, one can extract a correction function for the SiPM non-Unear 
response as the residual to linearity of the inverted SiPM response function, 

p,x; 1^^^^^ 

In the analysis presented here, we do not directly apply this formula but use tabulated mea- 
surements of A/pix versus Np^ for each SiPM instead. We fit the tabulated measurements with a 
double exponential function and use this to extract the value of /sat(A'pix) applied in the calibration 
chain (Eq. ^). The comparison of tabulated AHCAL data with the single and double exponential 
functions are shown in figure ^ for one example SiPM. 

The use of a second exponential in the fitting function does not have a solid physics motivation 
yet, tough it could be explained thinking of two areas of the SiPM active surface differently illu- 
minated from the WLS fiber. The total number of pixels is divided in two groups, and each of the 
group is described by an individual exponential function, such that the fitting function is the sum of 
two exponential functions like that in equation ^j^. There is no reson for the groups to be exactly 
two; this method could be extended to more exponentials. For practical reasons though, one needs 
to limit the number of free parameters. The choice shown in figure Q turned out to be sufficiently 
accurate and stable. The correction factor is close to unity for signals of about 30 pixels or 2 MIPs, 
and increases exponentially up to infinity for signals in saturation. 
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As illustrated for one SiPM in figure 0, the response curves of each SiPM has been sam- 
pled with 20 measurement points on a test bench setup illuminating each SiPM with LED light of 
variable intensity. For these studies, the SiPMs were not mounted on tiles, but were bare SiPMs. 
Therefore, all the pixels have been illuminated with hght in a homogeneous way. The measure- 
ment results for all SiPMs installed in the AHCAL are given in The maximum number of 
fired pixels (NtotQjaie)) for each SiPM is extracted with a fit to the measured points using Eq. 3.4. 
The spread (RMS) in the values of Ntot{bare) between all the curves is about 20 %. SiPMs with 
Niot{bare) > 900 have been pre-selected. This ensures not too large variations in the non-linear 
response function of each device. The 20 measurement points for each SiPMs are stored in a 
database. A linear interpolation of these points is used to calculate /sat from Eq. 3^ and linearize 
the calorimeter response during data reconstruction. 

Alternatively, Mot has also been extracted using the AHCAL LED monitoring system from 
measurements with the with SiPM mounted on a tile (A/^totCmounted)). As the saturation point 
in number of pixels is independent on the linearity of the light, no correction of the LED light 
intensity with PIN diode has been applied for this study. Figure ^ shows the ratio of A/^totCmounted) 
to MotCbare). The plot shows that the maximum number of pixels in the in-situ setup is on average 
80.5 % of the value determined in the laboratory setup |16] with a wide distribution (RMS=9 %). 
This factor is interpreted as geometric mismatch between the WLS fiber and the photodetector. 
The fiber has a 1mm diameter while the SiPM active surface area is Ixlmm^; the geometric 
ratio between areas is 79 %, in agreement with the measured value. Therefore, only a fraction of 
the SiPM surface is illuminated and the laboratory curves are re-scaled by the measured value of 
80.5 % to correct for this effect before they are used to correct for the SiPM saturation. 

The uncertainty of the determination of the saturation point for a single channel is lower than 
3 %, if the LED light range properly covers the SiPM saturation region, and if this region is mea- 
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Figure 6. Top view of the CERN beam test setup. The plot shows the instrumentation in 2007 (the y-axis 
is not to scale). The beam enters from the left side. See text for explanations of the components. 

sured well below the ADC saturation. Unfortunately, these conditions are true only for a sub- 
sample of about 73 % channels. Also judging from the tails in the distribution of Figure § some of 
the fitted results need to be investigated more accurately. For this reason, an average scaling factor 
is used for all channels. Further studies will address the possibility of using a channel-by-channel 
factor instead. Furthermore, the measured SiPM response points, from which the correction of non- 
linear detector response is calculated, are affected by the SiPM gain uncertainty of 2 %, discussed 
in the following section. 

4. The test beam experiment 

4.1 The experimental setup at CERN 

The data discussed in the following were collected in July 2007 at the CERN SPS test beam facility 
H6. A sketch of the experimental setup is shown in Figure |6[ Apart from the fully equipped 
AHCAL and a prototype of a tail-catcher and muon tracker (TCMT JI^]), the beam installation 
consists of various trigger and beam monitoring devices. A threshold Cerenkov counter was used 
to discriminate between electrons and pions. The beam trigger was defined by the coincidence 
signal of two plastic scintillator counters with 10 x 10 cm^ area, referred to as Scl and Sc2 in 
Figure |6[ One scintillator trigger (VI), with an area of 20 x 20 cm^ and analog read out, tagged 
multi-particle events. Another scintillator with a 100 x 100 cm^ surface and a 20 x 20 cm^ hole 
in the center (V2), was used to reject the beam halo. Three delay wire chambers (DCl, DC2 and 
DC3) were used to monitor the beam and reconstruct tracks. Events tagged by a scintillator with 
100 X 100 cm^ area (Mcl), placed behind the TCMT are most likely to be muons. 



During most of the tests, a silicon tungsten electromagnetic calorimeter [ ]19| ] was placed in 
front of the AHCAL, but this was not the case for the results reported here. The AHCAL was 
placed on a movable stage, which could shift the detector vertically and horizontally. In addition, 
the detector can be rotated with respect to the beam direction from an angle of 90° (beam normal 
to the detector plane) to approximately 60°. 

4.2 Monte Carlo simulation 

The test beam setup as shown in Figure |6| is simulated with Mokka [20], a GEANT4-based [13] 
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Monte Carlo program, followed by a digitization package simulating the response of the detector 
and electronics. The particle source of the simulation is positioned upstream of the Cerenkov 
detector. The beam position and spread are chosen to match the beam shapes measured in data by 
the delay wire chamber, DCS. The beam particles are parallel to the beam axis, according to the 
measurements in the three delay wire chamber detectors. The material upstream of the AHCAL is 
simulated. The sub-detectors are simulated with different levels of detail, depending on their impact 
on the physics analysis: material simulation only for the Cerenkov counter, raw energy depositions 
stored for the trigger counters, and partial electronics simulation for the tracking detectors. For the 
AHCAL, the simulation gives the raw energy depositions in a virtual scintillator grid of 1 x 1 cm^ 
tile size. The simulation is followed by a digitization procedure, which takes into account 

• the reaUstic detector granularity, 

• Ught cross-talk between neighboring tiles, 

• non-linearity and statistical fluctuations on the pixel scale, 

• SiPM and readout electronics noise. 

The actual geometry of the AHCAL is simulated by summing up the signal yield of 9 (36, 144) 
virtual cells to obtain those of the actual geometry 3 x 3 (6 x 6, 12 x 12) cm^ cells. 

Light cross-talk between neighboring cells, due to the imperfect reflective coating of the tile 
edges, is simulated assuming that from each 3 cm-long tile edge 2.5 % of the scintillator Ught leaks 
homogeneously to the neighboring tile. This value is scaled to take into account the fraction of edge 
shared with the neighbors for cells of different size. The amount of light cross-talk was checked 
experimentally only for two tiles. The leakage from one tile edge was quantified to be about 2.5 %. 
No information on the spread of this value between all tiles is given. This value is expected to 
influence the energy reconstructed and the transverse shower profile. From the comparison of the 
energy reconstructed in simulation and data, the value of 2.5 % for the light cross-talk on each tile 
edge is found to be adequate. A light cross-talk of 1.25 % or 3.75 % leads to a difference in the 
energy scale between data and Monte Carlo larger then 5 %. 

To simulate the non-linear behavior of the photodetectors, the energy deposition is translated 
from GeV to the number of fired SiPM pixels. For this, an intermediate step converts the response 
simulated in units of GeV to MDP equivalents. The conversion factor is estimated from the simu- 
lation of an 80 GeV muon beam in the AHCAL and is found to be 816keV/MIP, corresponding to 
the energy lost by a minimum ionizing particle in the scintillator. The amplitude in units of MIPs 
is then converted into pixels, using the measured light yield for each individual channel. With 
this scale the measured SiPM response curves from the test bench are used to simulate the SiPM 
non-linearity. Where not available, the curve of the next neighboring tile is applied. 

If A^pix is the ampUtude in pixels obtained this way and A^max is the saturation level of the 
individual channel, statistical effects are accounted for by generating a binomial random number 
with A'^max repetitions and a probability of Npa/N^sLx- The result is treated as the number of pixels 
firing for this specific event, and is translated back to the MIP scale with the channel-specific light 
yield. 
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At this stage, the Monte Carlo signal simulated the response of the AHCAL to the energy 
deposited by particles in an event. However, both the electronic components and the SiPM dark 
current induce noise. This noise component is assumed to be completely independent of the physics 
signal in each channel. The noise distribution is non-Gaussian due to the positive SiPM dark noise 
component. As the frequency of dark noise and the amount of inter-pixels crosstalk varies from 
SiPM to SiPM, the best way to account for noise is to take it from the data rather than try to 
simulate it. Noise events for each calorimeter cell are taken from data, are calibrated to the MIP 
scale according to the reconstruction procedure described in section ||, and are added randomly to 
the deposited energy of that given cell of a simulated event. 

A cell that could not be calibrated in the real detector, either due to an inactive photodetector 
or to missing calibration values, is also ignored in the simulation. This is about 2 % of the total 
number of cells in the calorimeter. 

5. Calorimeter response to positrons 
5.1 Selection of positron events 

The analysis presented here is based on positron runs between 10 and 50GeV. Each energy point 
has more than 150k recorded beam triggers. All positron runs have been simulated with statistics 
similar to the corresponding data runs. 

Single positron showers are selected for analysis using the beam instrumentation. Although the 
beam configurations are set to deliver a positron enriched beam, some contamination, mainly from 
muons, exists. The pion contamination is expected to be negligible, since the tertiary positron beam 
is produced from a higher-energy mixed beam impinging on a thin (2 Xq) lead target which does not 
result in the production of lower-energy tertiary pions. The muon contamination originates from 
in-flight decay of hadrons upstream of the production target, which results in a muon component 
that passes the momentum selection. 

Cells with a signal above threshold are called hits and E^n > 0.5 MIP is required. To reject 
empty events that can occur due to random triggers or scattered particles, the number of hits has to 
be A^hit > 65. Furthermore, the energy weighted center-of-gravity in the beam direction (z), defined 
as (z) = T.iZiEi/Y,iEi, has to be (z) < 390 mm (about half of the calorimeter depth). This require- 
ment eliminates muons, which deposit their energy equally distributed over the entire calorimeter 
depth, as opposed to electrons which have a short shower contained in the first half of the calorime- 
ter. It was found that this muon rejection was more efficient than the selection based either on 
the Cerenkov counter, which does not provide electron-muon separation for 30 GeV and above, 
or on the muon trigger Mcl which has an efficiency of about 50%. Particles which interact in the 
material upstream of the AHCAL aie removed by requiring a good track in the delay wire cham- 
bers (;^^/dof < 6), and a MlP-like energy deposition in the multiplicity counter (VI). With these 
selection criteria, 45 % of all recorded events at 10 GeV are accepted. According to Monte Carlo 
studies 99.9 % of all electron events pass the selection criteria, whereas 99.8 % of all muon events 
are rejected. The typical fraction of muons in a run is about 5-10 %. 

The uncertainty on the mean energy of the beam is reported in [ p| ] to be 
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Figure 7. Hit energy spectrum for 40 GeV positron showers compared to that of 40 GeV and 80 GeV pion 
showers from a Geant4 simulation. 

The first term is related to hysteresis in the bending magnets while the calibration and the un- 
certainties on the collimator geometry give the constant term. Since this uncertainty is negligible 
compared to the detector uncertainties, we assume the beam energy to be fixed. The dispersion of 



the beam energy can be calculated according to [ ]22| ] from the settings of the momentum selecting 



collimators on the beam line and is below 0.24 % for all the runs in this analysis. 
5.2 Linearity 

The linearity of the calorimeter response for a large range of incident particle energies is a key 
feature, which allows for an important test of the calibration chain. Electromagnetic showers offer 
the most rigorous test for non-linearity correction, since the energy deposited per single tile in 
an electromagnetic shower is larger than that in a hadronic shower for the same particle energy. 
Figure ^ shows the hit energy spectrum of a 40 GeV positron shower compared to the spectra of 
40 GeV and 80 GeV pion showers. The positron shower clearly has more hits with high energy 
deposition, even when the total particle energy is only half that of the pion. 

A set of positron runs with incidence normal to the center of the calorimeter face is analyzed. 
To minimize the influence of noise, the energy is summed up in a cylinder around the shower axis, 
where the shower axis is defined by projecting a track formed in the tracking system into the first 
layer of the AHCAL. This cylinder, sketched in Figure |8| (left), has a radius of 5 Moliere radii 
(r = 5 Rm, with Rm = 2.47 cm [j8|]), which ensures a lateral containment of more than 99 % of the 
shower energy. 

The length L of the cylinder is chosen to contain the whole shower energy. As suggest by 
simulations of 50 GeV electron showers, setting L to 20 layers contains the showers. Figure || 
(right) shows the final reconstructed spectra for positron runs in the energy range 10 to 50 GeV. 
The positrons are normally incident on the calorimeter front face, with a distribution centered in 
the same calorimeter cell for each run. The distribution is fit with a Gaussian function in the 
range ±2a. The position of the peak is taken as the mean energy response, ^mean. measured in 
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Figure 8. The shower energy is summed up in a cyhnder (left); see text for details. Spectra of the energy 
sum for positron data with energy between 10 GeV and 45 GeV (right). For each spectrum the mean energy 
response in units of MIP, Emean, is obtained with a Gaussian fit in the range ±2cT. 
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Figure 9. Reconstructed energy of a 1 GeV positrons for data (dots) and for Monte Carlo (filled histogram), 
as well as a Gaussian fit to data (blue line). 

units of MIP. The reconstructed energy of a 10 GeV positron shower is compared to the digitized 
energy from a Monte Carlo simulation in Figure ^ The agreement between data and simulation is 
satisfactory. 

The statistical uncertainties on the mean energy deposition are negligible. The main source of 
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systematic uncertainties is 2 % from MIP scale calibration. The uncertainty of 2 % on the SiPM 
gain determination, resulting from the fit stability and the uncertainty on the determination of the 
SiPM saturation level both affect the correction of the SiPM non-linear response. For the satura- 
tion level a common re-scaUng factor is applied to all SiPM curves determined in the laboratory 
setup. The rescaling is needed to account for the partial illumination of the SiPMs from the WLS 
fiber as discussed in Section 3.2. As shown in Figure ||, the ratio between the in-situ measured 
SiPM saturation level and the test-bench determined value has a wide distribution. Since a com- 
mon factor of 80.5 % is used to rescale all SiPM response curves, an uncertainty of 11.3 % on 
this value is assumed, which represents the spread of all measured values as taken from Figure 
To account for this uncertainty in an uncorrelated way for all SiPMs, 100 experiments have been 
performed assigning different rescaling coefficients for each channel, generated randomly with a 
Gaussian distribution centered at 0.80 and with a sigma of 0.09. For each experiment the energy in 
the calorimeter is reconstructed, using the set of curves rescaled by these randomly generated co- 
efficients to correct the non-linear SiPM response. Finally, the one standard deviation spread of the 
100 reconstructed energies from these simulated experiments is taken as the systematic uncertainty 
for the reconstructed energy. All of the above listed systematic uncertainties are uncorrelated and 
thus added in quadrature. The total systematic uncertainty ranges from 0.2 GeV (2 %) at 10 GeV 
to 1.7 GeV (3.4 %) at 50 GeV. 

The reconstructed energy in GeV is obtained as freco = £^mean/w^> where w is the electromag- 
netic energy scale factor (MIP-to-GeV). The scale factor is determined with a linear fit from zero to 
50 GeV to the distribution f'meaniMIP] versus f'beamlGeV]. The resulting values for data and Monte 
Carlo are Wdata = (42.3 ±0.4) MIP/GeV and wmc = (42.0 ±0.4) MIP/GeV, respectively. Within 
the uncertainties, the scale factors are in good agreement. 

The linearity of the AHC AL response to positrons is shown in Figure 10 A comparison of 
the data before and after correction for the SiPM non-linear response indicates the magnitude of 
this correction, which does not exceed 10 % even at 50 GeV positron energy. The values shown in 
Figure |l^ are reported in Table [I| 

The residuals for data and Monte Carlo are presented in Figure |ll|. Here, the green band 
indicates the quadratic sum of the energy dependent systematic uncertainties. In Table 1 the con- 
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Table 1. AHCAL energy reconstructed in data and MC (in units of GeV) for various positron beam energies. 
The table reports the values plotted in Figure 10. The systematic uncertainties for data are detailed in their 
percentage values. The total absolute error A|^' is the sum in quadrature of the uncertainties on the MIP, on 
the SiPM gain and on the saturation point determination. 
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Figure 10. Linearity of the AHCAL response to positrons in the range 10-50GeV. The blue dotted Une 
shows exact Hnearity. Dots correspond to data corrected for SiPM non-Hnear response, blue triangles show 
the data before this correction, and the open triangles show the simulation. The green band indicates the 
systematic uncertainty as quoted in Table 1, A^' [GeV]. 



tribution to the uncertainty from the SiPM gain variation, 5^^'", and from the saturation point 
determination, 51^" are hsted. The uncertainty on the MIP scale, 5^^^, cancels in the ratio since the 
same calibration constants are used in data and Monte Carlo. In Figure 11 (left), the residuals from 
the linear function suggests a non-zero offset at zero energy. This negative offset is the combined 
effect of the 0.5 MIP threshold (loss of energy) and the detector noise (addition of energy). Instead 
of the more conventional linear function with b = 0, the function ^mean = « • ^^beam + b can be used 
to fit the data in the range 10-50 GeV. A value of b = -10.3 ±7.4 Me V is found for the Monte 
Carlo offset. Once this offset is removed the Monte Carlo linearity is better than 0.5 % over the 



whole range, as shown in the right plot of Figure 1 1 



The deviation from linearity (Fig. 1 1 left) in data is less than 1 % in the range 10 to 30 GeV 
and the maximum deviation is about 3 % at 50 GeV. The remaining non-linearity at high energies 
hints at problems with the rescaling of the saturation curves, as described in Section 3.2. This 
behavior is not sufficiently reproduced in the Monte Carlo digitization, where the same curve is 
used to simulate saturation as is used to correct for it. 



The impact of the saturation correction is better seen in Figure 12 where the energy per hit 
is shown with and without the correction factor /sat applied, for 30 GeV electromagnetic showers. 
Whereas the correction is negligible for low signal amplitudes, it becomes significant at larger 
amplitudes, resulting in a strong correction for the tail of the distribution. The maximum energy 
deposited in one cell for a 30 GeV electromagnetic shower is ~230 MIPs corresponding to about 
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Figure 11. Residual to a fit of the data and Monte Carlo points presented in Figure using, the function, 
y = ax, (left), and the function, y ^ax + b (right), in the range 10-50 GeV. Dots correspond to data, and open 
triangles to simulation. The green band indicates the sum in quadrature of the energy dependent systematic 
uncertainties, 5^ and in Table 1 . 

3450 pixels (assuming a light yield, LY =15 pixel/MIP). For this amplitude the correction factor is 
/sat(^() ~ 3.1. The remaining miss-match between data and Monte Carlo around 100-200 MIPs is 
an effect of the non-perfect correction of the non-hnear SiPM response. This imperfect correction 
affects only a small fraction of the total energy; the hits above 50 MIPs contribute only 0.5 % (4 %) 
of the total energy at 10 GeV (40 GeV). 

5.3 Electromagnetic energy resolution 

Energy resolution is a principal figure of merit in calorimetry and is estimated as the width divided 
by the mean of a Gaussian fit to the energy sum within it 2 a of the mean of an initial fit over the 
full range. The resolution achieved with the AHCAL is plotted as a function of the beam energy 
in Figure 13. The values shown in this figure are reported in Table |2[ Fitting the AHCAL energy 



resolution in a range of ±2 a, with 

Oe a c 



results in a stochastic term of a = (2 1 .9 ± 1 .4) % \/E [GeV], whereas the constant term is = ( 1 .0 ± 
1 .0) %. The noise term of c = 58.0 MeV is extracted from the spread (RMS) of the random trigger 
event distribution and kept constant during the fitting procedure. The energy resolution agrees well 



with that of an earlier prototype (Minical) with 108 channels and of the same sampling [24], that 



was tested in the energy range between 1 and 6 GeV and reached a resolution with a stochastic term 
a = (20.7 ±0.7) %y^£'[GeV] and a constant term b = (2.6 ± 1.3) %. 

The energy resolution of the simulation is found to have a stochastic term of a = (21.5 it 
1 .4) % s/ E [GeV] , a constant term oib = (0.7 it 1 .5) % and again a fixed noise term of c = 58.0 MeV. 
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Figure 12. Hit energy spectrum for 30 GeV positron showers in the AHCAL. Open circles (dots) show the 
data before (after) correction for the non-Hnear response of the SiPM. The left plot shows the hit distribution 
in a logarithmic scale and the right plot on a linear scale. The shaded histogram is from digitized simulation. 

Within the fit uncertainty, the stochastic terms of data and simulation are in good agreement. The 
noise term is fixed to the same value as for data since the noise in the simulation is artificially added 
from random trigger data events. The constant term b, representing calibration uncertainties and 
non-linearities, is consistent with zero in the simulation as expected, since the same curves are both 
in the simulation of the non-linear SiPM response and in its correction. 
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Table 2. AHCAL energy resolution in data and MC for various positron beam energies. The table reports 
the values plotted in Figure 13. The listed uncertainties include statistical uncertainties and systematic 
uncertainties added in quadrature. 



5.4 Shower profiles 

The longitudinal profile of a shower induced by a particle with incident energy E in GeV traversing 
a matter depth t can be described as [ ^ ] 
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Figure 13. Energy resolution of the AHCAL for positrons (dots). The resolution agrees with that of a 
previous prototype (full triangles) with the same sampling structure. The errors are the quadratic sum of 
statistics and systematic uncertainties. The open triangles are the obtained from the analysis of the digitized 
simulated events. Fit curves to the data and MC are shown in the region 10-50 GeV. The dashed line is the 
extrapolation of the fit to AHCAL data in the low energy region covered by the MiniCal data. 



where the parameter a is an overall normalization, and the parameters CO and b are energy and 
material-dependent. The first term represents the fast shower rise, in which particle multiplication is 
ongoing, and the second term parametrizes the exponential shower decay. Given this parametriza- 
tion with t in units of radiation lengths, the particle multiplication and the energy deposition reach 
their maximum after 



0.5 



(5.4) 



radiation lengths from the beginning of the cascade of a particle with energy E. The critical energy, 
£c is a property of the calorimeter material and does not depend of the energy of the particle. The 
position fniax is called the shower maximum. 

The mean longitudinal profile of a 10 GeV positron shower is shown in the left plot of Fig- 
ure [I^ Due to the high longitudinal segmentation of the AHCAL, the shower rise, maximum and 
decay are clearly visible. Data and simulation are in qualitatively good agreement. To quantify 
this agreement, the profiles at each recorded beam energy are fitted with Eq. 53 and the maximum 
shower depth calculated as f^ax = co/b. The development of the shower maximum as a function of 
the beam energy is shown in the right plot of Figure 14 . The error bars show the uncertainty from 
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Figure 14. Longitudinal profile of a lOGeV positron shower in units of Xq (left) and scaling of the shower 
maximum as a function of the incident energy (right). The reconstructed energy (left plot) is shown for data 
(solid points), simulation (shaded area) and a fit to the data using Eq. 4.2 (line). The bottom insert shows the 
data/Monte Carlo comparison. The shower maximum (right plot) is shown for data (dots), simulation (open 
triangles) and the theory expectation given in Eq. 4.3 (solid line). 



the fits. The extracted shower maxima of both data and simulation are in good agreement with the 
theoretical behavior for a pure Fe calorimeter with a critical energy, from [23|, of £c = 21.04MeV, 
given in Eq. 5.4. 

The transverse shower profile of a 15 GeV positron shower is shown in Figure |l^ together 
with a simulation. The radius, p , is calculated with respect to the track of the incoming particle 
extrapolated from the tracking system to the AHCAL front face. Therefore, the radius is defined 
as pf = {xi — .Xtrack)^ + {ji " 3'track)^> whcrc {xi,yi) are the coordinates of the calorimeter cell with 
signal above threshold. The energy deposited in a calorimeter cell is normally assigned to the 
center of the cell. For the radial profile studies it is redistributed uniformly in bins of 1 mm^ 
before being assigned to one annular bin of inner radius p. In this way the energy deposited in one 
calorimeter cell can be shared between two adjacent annular bins. Proper normalization accounts 
for the fraction of the calorimeter cell area covered by each annular bin. The data indicate a 
broader shower than expected from simulation. The calculated mean shower radius ((/?) = 
for 15 GeV showers in data is about 9 % larger than the simulated one. 



The energy dependence of {R) is shown on the right plot of Figure 16 (left). The difference 
is almost energy independent. For completeness also the comparison of the RMS (y(^^)^^(^p) 
of the shower radius distribution is shown in Figure [T^ (right). An energy dependent disagreement 
of data and Monte Carlo is observed for this variable which increases to a maximum of 7 % for 
50 GeV. Several studies have been performed to find the cause of this effect including noisy and 
inactive cells, different beam shape, influence of the light cross-talk between tiles, misalignment 
between calorimeter and tracking system and of calorimeter layers. The broader shower in data 
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Figure 15. Transverse profile of a 15 GeV positron shower. The energy density is shown in 10 mm wide 
concentric rings centered around the shower axis. 



is still not understood and further studies of asymmetric light collection on the tile, the influence 
of varying dead space between tiles due to varying thickness of reflector coating, etc., will follow 
to investigate the discrepancy. For the purpose of the validation of the calibration procedure the 
current level of agreement is acceptable, though this mismatch will have to be taken into account 
when comparing hadronic shower shapes. Furthermore, hadronic showers have a much smaller 
energy density than electromagnetic showers; therefore, any local effects, (i.e. the impact of dead 
areas or misalignment between layers), are strongly amplified in electromagnetic showers, while 
the influence is expected to be much less pronounced in hadronic showers. 
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Figure 16. Mean (left) and RMS (right) of the transverse shower distribution as a function of beam energy. 
Dots are from data and open triangles are from simulation. 
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Figure 17. Schematic view of tile positions in an AHCAL scintillator plane used for the uniformity test 
(left) and uniformity of the calorimeter response for various positions of incident beam with respect to the 
detector (right). Tile position eight is approximately in the center of each calorimeter layer. The dashed lines 
show the systematic uncertainties. Statistical uncertainties are negligible. 



6. Uniformity studies 

6.1 Uniformity of the calorimeter response 

The uniformity of the AHCAL response is explored by shifting the AHCAL to different positions 
with respect to the beam axis, at normal angle of incidence. This procedure is visualized in the left 
part of Figure 17. Each square in the sketch represents one scintillating tile of 3 x 3 cm^ and beam 
events with a track pointing to a 1 x 1 cm^ region centered on each tile in turn were selected. The 
uniformity of the calorimeter response at the 15 different positions has been tested. For this study 
10 GeV positron runs are analyzed, where the movable stage was used to displace the calorimeter 
in the x-y position with respect to the beam-line (z-axis). 



As shown in the right plot of Figure [17|, when excluding position 10, the uniformity of the 
calorimeter response is better than 2.1 %. The 10% deviation between reconstructed and beam 
energy at position 10 is due to an inactive cell at the shower maximum, which is not corrected in 
the calibration. 

6.2 Angular dependence of the calorimeter response 

The movable stage carrying the AHCAL was used to collect positron data at incident angles of 
90°, 80°, 70° and 60°. The rotation and staggering of the AHCAL are sketched in the left plot of 



Figure 18, where the beam is entering from the top. In the rotated configuration, the modules were 
staggered to ensure the highly granular core of 3 x 3 cm^ was aligned with, and hence sampled, the 
shower core. 

For each angle of incidence 6 several 10 GeV positron runs with different impact points on 
the calorimeter front surface were taken. The average of all the runs is used to define Eyeco at one 
given 6, while {E,-eco) is the mean reconstructed energy at any angle of incidence. The normalized 
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Figure 18. Schematic view of the AHCAL rotated with respect to the beam (left) and reconstructed energy 
of 10 GeV positrons normalized to the average versus angle of incidence (right). To improve legibility, the 
data (solid points) and the simulation (red triangles) are slightly shifted in opposite directions on the abscissa. 
The systematic uncertainty is shown by dash-dotted hnes. Additionally, the spread of all measurements 
performed at one inclination angle are shown as an error for each point. 



reconstructed energy is plotted in Figure |18| as a function of 6. The spread (RMS) between the 
various analyzed runs per inclination angle is used as the systematic uncertainty. This spread 
is smaller than the calibration systematic uncertainty in the calorimeter, shown in the plot as an 
error band around the ratio of one. Showers at various inclination angles only partially share the 
same calorimeter cells, therefore the full systematic uncertainty from calibration is an overestimate 
of the real error, but the spread between measurements performed at one inclination angle is an 
underestimate. Taking this into account, the increase in the reconstructed energy of data between 
60° and 90° is not significant. A more precise analysis would require more data at different angles 
which are not available at present. 

6.3 Influence of cell structure 

The scintillating tiles used in the AHCAL have a WLS fiber embedded in a groove, a SiPM inserted 
into a small groove on one end of this fiber, and a mirror in a groove on the other end. This structure 
is visible in the picture of a 3 x 3 cm^ tile in Figure 19, where the SiPM is located in the lower left- 
corner, the mirror is placed in the diagonal corner, and the WLS fiber is embedded in a quarter 
circle. 

Since electromagnetic showers have a short transverse extension, the impact of the cell struc- 
ture slightly reduces the energy resolution of the calorimeter. To study this effect, we take advan- 
tage of the delay wire chambers that were present in the beam line. They are used to reconstruct 
the track of the incoming particle. This track is then extrapolated to the front face of the AHCAL. 
The shower energy (energy summed over the entire calorimeter) for 10 GeV positrons with this 
impact position, normalized to the shower energy averaged over all impact positions, is plotted in 
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Figure 19. Picture of the scintillating tile (left). Effect of the AHCAL scintillator tile structure on the energy 
measurement (summed over the entire calorimeter) for 10 GeV electromagnetic showers (right). 



Figure |T^. 

As shown in the figure, the measured energy drops slightly over the area of the WLS fiber. A 
particle with a trajectory intersecting the SiPM (in the lower left corner of the plot) or the reflecting 
mirror at the end of the WLS fiber (in the upper right corner of the plot) shows a significant loss of 
response with respect to the tile average by about 8 % and 4 % respectivley. At the position of the 
WLS fiber the tile response is about 2 % lower than average. The drop at the other two corners of 
the tile in this study reflects the energy loss associated with the SiPMs located in the neighboring 
tiles as the observable used is the energy summed over the entire calorimeter. Measurements of 



single tile uniformity using a collimated source have been performed and are reported in [ 25 



These measurements confirm a lower response of electromagnetic showers hitting the SiPMs or 
the reflecting mirrors. Though this large degradation (8 % at the locations of the SiPMs) is quite 
unrealistic in a collider detector, where the particles are always traversing the calorimeter under an 
angle. In this case the tile response non-uniformity averages out with no influence on the energy 
resolution. Furthermore, electromagnetic showers have a short lateral extension. For pion showers, 
which are much wider, the effect has not been observed in data. The effects of gaps between the 
calorimeter tiles, as well as the non-uniform response of the tiles, in view of the impact on the 
energy resolution, have been studied using Monte Carlo events. The results are reported in [27] 
and show that these type of effects do not have a significant influence on the measurement of hadron 
showers. 



7. Conclusions 

The response of the CALICE analog hadron calorimeter to positrons was measured for energies 
between 10 and 50 GeV, using data recorded at CERN in summer 2007. The calorimeter response 
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is linear to better than 3 %. A better SiPM saturation correction would improve the linearity, and 
for future developments a larger dynamic range is desirable. This study is ongoing, but the ef- 
fect on pion energy reconstruction will be negligible due to the much smaller energy per hit in a 
hadronic shower compared to an electromagnetic shower. The energy resolution for positrons is 
found to have a stochastic term of (21.9 it 1.4) %\JE [Ge V] , and a constant term of about 1 %. 
Good agreement between data and simulation validates the simulation of the various detector char- 
acteristics. For comparison, Ref. [28] reports for the ATLAS tile calorimeter an energy resolution 
of 28 %^JE [GeV] stochastic and 2.8 % constant term for electrons at 20 deg from normal inci- 
dence. This is also and hadron sampling calorimeter alternating steel and scintillator tiles, but with 
a much coarser granularity than the AHCAL and a different a readout via standard photomulti- 
plier tubes. The same readout technology as in the AHCAL is also implemented in a scintillator- 
Tungsten electromagnetic calorimeter, ScECAL build within the C ALICE collaboration [29]. For 
this calorimeter the energy resolution to electrons is of (15.15 ±0.03) %\JE [GeV] stochastic and 
(1.44 lb 0.02) % constant term. This analysis provided confidence that the detector performance 
and simulation are sufficiently understood to pursue the investigation of hadronic showers. 

Systematic studies are performed to investigate the quality of the calibration in as many 
calorimeter cells as possible. The uniformity of the calorimeter response to electromagnetic show- 
ers is studied with beams at different impact points and different incident angles. The results are 
consistent with no angular and spatial dependence within the quoted systematic uncertainty on the 
calibration procedure. 

The high segmentation of the AHCAL is well-suited for studying the longitudinal shower 
development with high accuracy and for determining the shower maximum. The point of maximum 
energy deposition along the shower propagation axis is located between 5.5 and IXq for the 
range of particle energies used, consistent with simulation and theoretical prediction. 

The transverse shower spread is more difficult to measure because it is sti^ongly affected by 
uncertainties in the beam profile, in the variation of light cross-talk between tiles, and in the mis- 
alignment of calorimeter layers. Currently, the data indicate a broader shower than expected from 
simulation. However, the level of agreement is acceptable for the validation of the calibration pro- 
cedure if one considers that the effect on hadronic showers will be less important due to the lower 
energy density of hadronic showers. 
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